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Abstract. A classical spin network consists of a ribbon graph (i.e., an abstract graph with a cyclic ordering 

_ of the vertices around each edge) and an admissible coloring of its edges by natural numbers. The standard 

^_>^ ' evaluation of a spin network is an integer number. In a previous paper, we proved an existence theorem for 

(N ■ the asymptotics of the standard evaluation of an arbitrary classical spin network when the coloring of its 

^ , edges are scaled by a large natural number. In the present paper, we extend the results to the case of an 

^ ■ evaluation of quantum spin networks of arbitrary valency at a fixed root of unity. As in the classical case, 

' our proofs use the theory of G-functions of Andre, together with some new results concerning holonomic 

' and g-holonomic sequences of Wilf-Zeilberger. 

00 
(N 

H 
O 



Contents 



1. Introduction 2 

1.1. History 2 

1.2. Classical spin networks 2 

1.3. Quantum spin networks 2 
CO' 1.4. Sequences of Nilsson type and G-functions 3 
^ \ 1.5. Statement of our results 4 

. I Q Acknowledgment 5 

, 2. Evaluation of quantum spin networks 5 

■ 2.1. Local projectors 5 
2.2. The standard evaluation of a quantum spin network 5 

' 2.3. Evaluation of spin networks using the shadow formula 7 

O ! 2.4. Integrality 9 

3. Proof of Theorems 1.3, 1.4, 1.5 and 1.6 10 
3.1. Behavior of q-holonomic sequences under differentiation 10 

, , ■ 3.2. Evaluations of q-holonomic sequences at roots of unity are holonomic 11 

/S[ 3.3. Proof of Theorem 1.6 13 

■ 3.4. Proof of Theorem 1.3 14 

4. Examples 14 

4.1. Evaluations of the q- multinomial coefficients at a root of unity 14 

4.2. The q-difference equation of the regular 6j-symbol 15 

4.3. The regular 6j-symbol at roots of unity 16 

5. Open problems 17 
References 17 



Date: November 3, 2010 . 

S.G. was supported in part by NSF. 

1991 Mathematics Classification, Primary 57N10. Secondary 57M25. 

Key words and phrases: Spin networks, quantum spin networks, ribbon graphs, 6j -symbols, Racah coefficients, angular 
momentum, asymptotics, G-functions, Kauffman bracket, Jones polynomial, Wilf-Zeilberger theory, holonomic sequences, 
q-holonomic sequences. 



1 



2 



STAVROS GAROUFALIDIS, STAVROS GAROUFALIDIS, AND ROLAND VAN DER VEEN 



1. Introduction 

1.1. History. Spin networks originally arose from calculations of angular momentum in quantum mechanics 
[VMK88]. They were formalized in the sixties by Penrose [PenTla, PenTlb] in attempt to quantize gravity 
combinatorially. Similar ideas were developed by Ponzano and Regge who related the semi-classical expansion 
of 3D gravity to the Regge action, [BL81a, BLSlb, Pen68, Wig41]. In recent years, spin networks have 
played an important role in the development of loop quantum gravity; see [BC98, EPR08, Per03, Rov04] 
and references therein. In the eighties, quantum spin networks were used by Reshetikhin, Turaev and Viro 
as building blocks of combinatorially defined invariants of knotted 3-dimensional objects; see [Tur94, TV92]. 
Those topological invariants that are collectively called TQFT generalize the famous Jones polynomial of a 
knot [Jon87]. 

A key problem in classical and quantum spin networks is their asymptotic behavior for large spins. In 
the case of the 6j-symbols Ponzano- Regge conjectured the leading term of an explicit asymptotic expansion 
and gave ample numerical evidence; see [Pen68]. The Ponzano- Regge conjecture was proven by Roberts 
[Rob02] using methods of geometric quantization. The existence of a general asymptotic expansion for all 
classical spin networks (to all orders in perturbation theory, in a constructive way) was recently obtained by 
the authors in [GvdV09] . The method of [GvdV09] was to convert the problem of asymptotic expansions to 
questions in Algebraic Geometry and Number Theory, and use the highly developed theory of G-functions. 
An important part in this conversion was the use of holonomic functions, developed by Zeilberger. 

1.2. Classical spin networks. A classical spin network consists of a ribbon graph F (i.e., an abstract 
graph with a cyclic ordering of the vertices around each edge) and an admissible coloring 7 of its edges by 
natural numbers. The standard evaluation (r,7) of a spin network is an integer number. In a previous 
paper [GvdV09], we proved an existence theorem for the asymptotics of the sequence (P, 717) of an arbitrary 
trivalent classical spin network when the coloring of its edges is scaled by a large natural number n. In 
addition, we presented several ways (algebro-geometric, combinatorial, and numerical) for computing the 
asymptotics of (F, nj) for large n. 

The goal of the present paper is to extend the results of [GvdV09] to the case of quantum spin networks 
at a fixed root of unity. As in the case of classical spin networks, our proofs use the theory of G-functions 
of Andre. There are two new ingredients which allow us to use the theory of G-functions. They come from 
the theory of holonomic and g-holonomic functions of Wilf- Zeilberger: 

(a) Theorem 1.4 which states that the derivative of a 5- holonomic sequence with respect to q is q- 
honolomic. 

(b) Theorem 1.5 which states that the evaluation of a g- holonomic sequence at a fixed root of unity is 
holonomic. 

1.3. Quantum spin networks. Quantum spin networks differ from their classical versions in two ways: 

(a) the underlying graphs F are knotted (i.e., embedded in S^) and not abstract, 

(b) their evaluations are polynomials of q, and not simply integer numbers. 

Recall that a knotted ribbon graph is an embedded framed graph in S^, of arbitrary valency, together with 
a cyclic ordering of the edges around each vertex. We will restrict ourselves to integer framings so that 
thickening the graph gives rise to an orientable surface. A quantum spin network consists of a knotted 
ribbon graph F together with a pair of functions 7 = (7£:,7y). Here je '■ Edges (F) — > N is an admissible 
coloring 7 of its edges, and 7y : Vert(F) — > {projectors} a choice of a local projector at each vertex of F. 
Local projectors are explained in detail in Section 2, and can be ignored in case F is a trivalent graph. 

The standard evaluation {T,^){q) of a quantum spin network is a rational function of q^^^', see [Tur94]. 
Recently, Costantino [Cos09] obtained an integrality result for the standard evaluation of a quantum spin 
network, namely: 

(1) (F,7)(<z)e<75z[(7±i] 

For some n depending on the network. Our first result deals with fixing the fractional power of q. 
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Proposition 1.1. For every quantum spin network (r,7) there exist a Z/AZi-valued quadratic form (3(7) 
such that 

(2) (r,7)(<7) = g'«(''\(r,7))(g), ((r,7)>(g) e nq^']. 

The above proposition allows us to define a modified evaluation ((r,7))(C) G C that can be evaluated a 
fixed root of unity ^. To make ((r,7)) well defined we use the convention that < Q(7) < 3. When C = 1, 
it follows from the definitions that ((r,7))(l) equals the classical evaluation and is hence independent of the 
embedding of T in 3-space. The goal of this paper is to extend the results of [GvdV09] to the case of the 
standard evaluation of a quantum spin network at a fixed root of unity 

1.4. Sequences of Nilsson type and G-functions. To state our results, we need to recall what is a G- 
function and what is a sequence of Nilsson type. Sequences of Nilsson type are discussed in detail in [GarlO]. 
We recall some definitions here for the convenience of the reader. 

Definition 1.2. We say that a sequence (a„) is of Nilsson type if it has an asymptotic expansion of the form 

(3) - -^"""(l0g«/^A,a,/3/iA,a,/3 

\,a.,l3 

where 

(a) the summation is over a finite set of triples (A, a, (3) 

(b) the growth rates A are algebraic numbers of equal absolute value, 

(c) the exponents a are rational and the nilpotency exponents /? are natural numbers, 

(d) the Stokes constants S\^a,/3 are complex numbers, 

(e) the formal power series h\^a,p{x) g ^[[a:]] are Gevrey-1 (i.e., the coefficient of is bounded by 
C"n! for some C > 0), 

(f) K is a number field generated by the coefficients of h\_ajj{x) for all A, a, (3. 

Definition 1.3. We say that series G{z) = X^riLo ^ G-function if 

(a) the coefficients a„ are algebraic numbers, 

(b) there exists a constant C > so that for every ri S N the absolute value of every Galois conjugate 
of a„ is less than or equal to C", 

(c) the common denominator of ao, . . . , is less than or equal to C", (where the common denominator 
c? of ao, . . . , On is the least natural number such that doi is an algebraic integer for t = 1, . . . , n), 

(d) G{z) is holonomic, i.e., it satisfies a linear differential equation with coefficients polynomials in z. 

The following connection between sequences of Nilsson type and G-functions was observed in [Gar09, 
Prop.2.5]. 

Theorem 1.1. [Gar09, Prop.2.5] If G{z) = X^i^o'^"'^" G-function, then (a„) is a sequence of Nilsson 
type. 

The above existence theorem has a valuable, effective corollary. 

Corollary 1.4. If G(z) = X^^o'*"-^" ^ G-function, and suppose that we are given either (a) a linear 
differential equation for G(z) with coefficients in L[x\ for some number field L, or (b) a linear recursion for 
(a„) with coefficients in L[n\, then one can effectively compute A, a, /3 and h\^a.p{x) G ^^^[[2^]] for a number 
field K such that the asymptotic expansion (3) holds, for some unknown Stokes constants Sx.ap- 

In other words, a linear recursion for (a„) or a linear differential equation for G(z) allows us to compute 
the asymptotic expansion (3) to all orders in 1/n, up to a finite set of unknown Stokes constants. For explicit 
illustrations of the above corollary, see [FS09, WZ85] and also [GvdV09, Sec. 10] where the authors discuss 
in detail effective asymptotic expansions of evaluations of classical spin networks. 
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1.5. Statement of our results. Let /''^•'(g) = d^'/dq^f{q) denote the r-th derivative of a Laurent polyno- 
mial fiq) e Z[g±i]. 

Theorem 1.2. For every quantum spin network (r,7), every complex root of unity ^ and every natural 
number r > 0, the sequence {{T,n"/))^^\C) of Nilsson type. 

Theorem L2 follows from Theorem 1.1 and the following theorem. 

Theorem 1.3. For every quantum spin network (r,7), every complex root of unity C, and every natural 
number r, the generating function 

oo 

(4) Fr,,,c,.(^) = E«r,"7))^'^HC)^" 
is a G-function. 

Theorem 1.3 follows from Theorems 1.4, 1.5 and 1.6 below, which involve structural properties of the 
classes of holonomic and g-holonomic sequences and may be of independent interest. To state them, recall 
that a sequence (/„) of complex numbers is holonomic if it satisfies a linear recursion of the form 

(5) Cd{n)fn+d H h ca(n)fn = 

for all n where Cj{n) G K[n\ are polynomials in n with coefficients in a number field K for j = 0, . . . , d with 
Cd 9^ 0. Likewise, a sequence {fnW)) of rational functions of q is q-holonomic if it satisfies a linear recursion 
of the form 

(6) Cd(g", q)fn+d{q) + ■■■ + c,{q", q)fUq) = 

for all n where Cj{u, v) G K[u, v] are polynomials in two variables u and v for j = 0, . . . , d with Cd ^ 0. 
Holonomic and q-holonomic sequences were studied in detail by Zeilberger; see [Zei90, WZ92]. 

Theorem 1.4. The derivative with respect to q of a q-holonomic sequence is q-holonomic. 

Theorem 1.5. For every q-holonomic sequence fn{q) £ ^[9^^] o,n-d every complex root of unity C, the 
evaluation /„(C) is a holonomic sequence. 

Let us make some remarks. 

Remark 1.5. Theorem 1.5 fails when C, is not a complex root of unity. For example, fn{q) = 9"^ is q- 
holonomic and satisfies the recursion fn+i{q) — Q^^~^^ fn{q) = 0. On the other hand, /n(a;) is holonomic only 
when w is a complex root of unity. Theorem 1.5 is another manifestation of the importance of roots of unity 
in Quantum Topology. See also Section 3. 

Remark 1.6. Theorems 1.4 and 1.5 presumably hold for for multi- variable g-holonomic sequences fni,...,nr (?) G 

Multi- variable holonomic and g-holonomic sequences were introduced and studied in [Zei90, WZ92]. 
In the present paper, we will not use them. However, the reader should keep in mind that for every quantum 
spin network (r,7) the multi- variable sequence 7 H> (T,^){q) is g-holonomic. This follows from Section 2.3 
where it is shown that (r,7)(g) is a multi-sum of a g-proper hypergeometric term. By the fundamental 
theorem of WZ-theory (see [WZ92]), (r,^){q) is g-holonomic in all 7-variables. 

Remark 1.7. The g-holonomic sequence /n(C) in Theorem 1.5 need not be exponentially bounded. For 
example, 

" 1 A: 

Uq) = n e ni] 

k=l ^ 

satisfies the linear recursion 

(g-l)/„+i(<7)-(9"+i-l)/„(9) = 0. 
Thus, fnio) is g-holonomic. Its evaulation at C = 1 is given by /n(l) = n\ which is holonomic, since it 
satisfies the linear recursion 

/„+i(l)-(?i + l)/„(l) = 
However, n! is not exponentially bounded. 
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Despite the above remark, the evaluations of quantum spin networks at a fixed root of unity is exponen- 
tiaUy bounded- 
Theorem 1.6. For every quantum spin network (r,7), every complex root of unity ^ and every natural 
number r, there exists C > (which depends on r,7 and r) such that 

(7) |((r,n7))('^)(C)l<C" 
for all n e N. 

1.6. Acknowledgment. The paper was conceived during a conference on Interactions Between Hyperbolic 
Geometry, Quantum Topology and Number Theory held in Columbia University, New York, in June 2009. 
The authors wish to thank the organizers, A. Champanerkar, O. Dasbach, E. Kalfagianni, I. Kofman, W. 
Neumann and N. Stoltzfus for their hospitality. In addition, we wish to thank F. Costantino, G. Kuperberg 
and D. Zagier for enlightening conversations. 

2. Evaluation of quantum spin networks 

2.1. Local projectors. We start by defining classical spin networks of arbitrary valency using the notion 
of a local projector. A classical spin network will be a pair (F, 7) of an abstract ribbon graph F (of arbitrary 
valency), together with a pair of functions 7 = (7B, 7y) where 7^ : Edges(F) — > N is an admissible coloring 
7 of its edges, and 7b : Vert(F) — > {projectors} a choice of a local projector at each vertex of F. A local 
projector is defined as follows. Given an admissible coloring (ci, . . . ,Cd) of the d edges around a vertex v 
of a ribbon graph, place a collection of ci + • • • + points on a disk. A local projector is a planar way to 
connect these points with disjoint arcs on the disk and no U -turns, as in the following example: 



Figure 1. An example of a projector for a degree 4 vertex with edge coloring (2, 5, 2, 3). 

Note that if p is a local projector and n is a natural number, then there is there is a canonical choice 
of a local projector np by cabling each arc oi p into n arcs. In [Kup96, p. 118], Kuperberg uses the term 
clasped web space W{ci, . . . , Cd) which has a basis the local projectors defined above. Local projectors are a 
pictorial way to encode dual canonical bases for 5(2, as was shown in [FK97]. However, clasped web spaces 
and dual canonical bases do not coincide for sts; see [KK99]. Finally note that there are alternative methods 
of defining multi-valent vertices, for example [Yok96] , [Bar03] . However our definition using projectors works 
without change for Lie algebras of arbitrary rank. 

2.2. The standard evaluation of a quantum spin network. Recall from the introduction that a quan- 
tum spin network is a pair (F,7) where F is a framed ribbon graph embedded in and 7 = (7£;,7v') as 
in the classical case. We allow F to have multiple edges and loops and it may be disconnected. Knots and 
links are quantum spin networks with zero vertices. However we restrict ourselves to integral framings of F. 
In this section we define the evaluation of quantum spin networks in terms of the Kauffman bracket. Recall 
that the quantum integer [n] and the balanced quantum factorial [n] ! is defined by 

A2n _ A-2n ^ 

(8) w - ^2 _ ^-2 ^ w!-nw 

where 

(9) = q. 



k=l 



6 



STAVROS GAROUFALIDIS, STAVROS GAROUFALIDIS, AND ROLAND VAN DER VEEN 



The Kauffman bracket evaluation is determined by the fohowing rules: 




QUD = -[2]-D 

To emphasize the similarity to the evaluation of classical spin networks we use the explicit expression for the 
Jones- Wenzl idempotent from [KL94, Chpt.3]. This is defined as the following formal sum of braids: 

Here b is the label of the edge, and for any permutation a we denote by the unique negative (with respect 
to orientation downwards) permutation braid corresponding to a. By i{a) we mean the minimal length of 
a written as a product of transpositions. Note that we leave out the quantum factorial [6]! that is used in 
[KL94]. 

Definition 2.1. (a) We say a quantum spin network is admissible if for every vertex v the projector 7y('y) 
matches the labels of the edges given by 7^;. 

(b) The evaluation (F, 7)^ of a quantum spin network (F, 7) is defined to be zero if it is not admissible. An 
admissible quantum spin network is evaluated by the following algorithm. 

• Use the ribbon structure to thicken the vertices into disks and the edges into untwisted bands. 

• Replace each vertex v by the pattern of the projector 7v(w) and replace each edge by the linear 
combination of braids as shown in Figure 2. 

• Finally the resulting linear combination of links is evaluated calculating the Kauffman bracket. 




Figure 2. The rules for evaluating a quantum spin network. Replace vertices and edges to get a linear 
combination of links that is evaluated using the Kauffman bracket. 

In the case A = —1 this definition agrees with the Penrose evaluation defined in [GvdV09]. For general A 
the definition agrees with that given in [MV94, KL94, CFS95] except for the missing quantum factorial [6]! 
in the denominator, one for every edge. The reason for leaving out this factor is that there is a better way 
to normalize (see Definition 2.2 below) which results in Laurent polynomials in q^^^ (as opposed to rational 
functions in q^^^) while the value of {T,n'y){q) at a fixed root of unity q grows at most exponentially with 
respect to n, see Theorem 1.6. 

Definition 2.2. The standard evaluation of a quantum spin network is defined by 

(r,7) = |^(r,7)^ 

Here [I]! = rii;[7v(w)]!i where for a projector [p\l means grouping all strands connecting the same two edges 
and forming the product of the quantum factorials of these numbers. 
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Note that with this normahzation subdividing an edge colored a into two edges eolored a connected by a 
vertex whose third edge is colored does not change the value of the evaluation. This is the way in which 
the value of the unknot should be interpreted in order to obtain the correct value (— l)°[a + 1] instead of its 
quantum factorial. 

We end this section by proving that any quantum spin network evaluation actually reduces to a trivalent 
quantum spin network evaluation. Hence all results previously known in the trivalent case extend to the 
general case. 

Lemma 2.3. Let (r,7) be a quantum spin network. There exists a trivalent quantum spin network (r',7') 
such that for all n e N 

(r,n7) ^ (r',n7') 

Proof. It is convenient to take a slightly different view of the evaluation algorithm defined in 2.1. Instead of 
expanding out the linear combination at every edge, we view the Jones- Wenzl idempotent as box with arcs 
coming out. After expanding the vertices the whole network becomes a number of boxes connected by arcs. 
The idea is to add extraneous boxes and to reinterpret the result as the evaluation of a trivalent quantum spin 
network. The key property that makes this work is the fundamental fact that once one includes the factor 
in a box with b strands, it becomes an idempotent [MV94]. Since we're not using this normalization we 
get that adding an extra box to b parallel adjacent arcs coming out of a box multiplies the evaluation by [b]l. 
This factor will cancel with the the normalization factor due to the trivalent vertices that will be created in 
the process. 

Let us first describe how to modify a single half-edge pointing into a multivalent vertex v with projector 
p. The arcs coming into v from half edge e are grouped into aj parallel arcs according to the projector 
p. To the box that already marks the beginning of e we now add an extra stack of boxes as shown in the 
middle of Figure 3. This will multiply the Penrose evaluation of F by a factor Oe = ([a2]!...[a„]!)([ai -|- 02... -I- 
a„_i]!...[ai 4- a2]![ai]!) 



ai 02 03 




Figure 3. Turning an edge into a trivalent network. 

Doing this for every half edge gives just enough boxes to split into a trivalent quantum spin network F' 
as shown in the above figure. However the smallest boxes at the tip of the tree meet another equal sized 
box on the other side which are not there in the evaluation of F'. Therefore we have to remove those boxes 
again, which reduces the Penrose evaluation by exactly a factor = [I{T)]\. Therefore, 



And hence passing to the standard normalization we find: 



(r,7> = 



[I]! 



= (r',7'> 



concluding the proof. 



□ 



2.3. Evaluation of spin networks using the shadow formula. In this subsection we describe a way of 
evaluating spin networks in terms of the shadow formula. We restrict ourselves to trivalent F. Since F is 
supposed to be orientablc, we can choose a blackboard framed diagram D of F. 

The shadow formula expresses the evaluation in terms of a multi-dimensional sum of Ij, 3j and 6j-symbols. 
The latter are the evaluations of three basic spin networks, shown in Figure 4. 
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Let 

(10) 



ai,a2, • ■ 



[ai]! 



denote the quantum multinomial coefficient when ai + . . . ~ a. The value of the Ij, Sj and 6j-symbols is 
given by the following lemma of [KL94] (see also, [MV94, Wes98]), using our normalization. 

Lemma 2.4. (a) We have 

{Q,a)^{^ir[a + 1] 

(b) Let (0,7) denote the 8 spin network admissibly colored by 7 = (a,6, c) as in Figure 4. Then we have 



(11) 



(6,7) = (-1) 2 [ 



2 

-a+b+c a—b+c a+b—c 
2 ' 2 ' 2 



(c) Let ( ,,:';> , 7) denote a tetrahedron labeled and oriented as in Figure 4 and admissibly colored by 7 
(a, 6, c, d, e, /). Then we have 

min Sj 

(12) (A' 7)= E (-i)'[fc + i] 



k 



Si — k, S2 — k, 5*3 — k,k — Ti, k — T2, k — T3, k — T4 

a; — iiicijt ±i 

where Si are the half sums of the labels in the three quadrangular curves in the tetrahedron and Tj are the 
half sums of the thee edges emanating from a given vertex. In other words, the Si and Tj are given by 



(13) 
(14) 



d + b + c) 



Ti 



1 



(a + 6 + e) 



To 



1 



(a + c + /) 



-d + e + f) 
1 



S^ = \{b + c 



{c + d + e) 



e + f) 

d + f). 



1 

2'" ■ " ' 2''' ' " ' 2'^ ' " ' 2 

In addition to the three basic spin networks above, we choose to consider a variant of the tetrahedron 
that represents a crossing. The crossing can be reduced to a 6j-symbol using the half twist formula, [MV94] . 



q 



ri(ri+2)-r3(T-3+2) + r2(r2+2) + r4(r4+2) 



We are now ready to explain how to evaluate a spin network using the shadow formula. A blackboard 
framed diagram of F gives rise to a planar graph D, whose edges arc colored by 7 and whose vertices are 
either vertices of F or crossings. Let V, E, F, C denote the sets of trivalent vertices, edges, faces and crossings 
(4-valent vertices) of D. We will express the evaluation of the quantum spin network as a sum over all 
admissible colorings of F by natural numbers. Here admissible means that for any two faces colored ri, e, r2 
and separated by edge e, the numbers (ri,r2,7(e)) satisfy the triangle inequalities. Moreover the outer face 
should get color 0. 



(r,7)=EnonQ-^nAn 



(15) 

r f£F e£E v^V c£C 

The sum is over admissible colorings r of the faces F and the labels of the symbols in the formula are found 
from the coloring and 7 as follows. The Ij symbol gets the color of the face it corresponds to. The theta 
gets the color of the edge it corresponds to plus the two colors of the faces it bounds. For every vertex the 
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Figure 5. An example spin network r,7 (left) and its diagram with encircled crossings and vertices and a 
coloring of the regions. 

corresponding 6j symbol is colored by the three adjacent faces and the three adjacent edge labels. Finally 
every crossing is counted by a skew 6j-synibol obtained by encircling the crossing and coloring according to 
the faces one crosses. For a more elaborate discussion see [Cos09]. 

Consider the following simple example of a spin network evaluation using the shadow formula. 

Note that one can find the contributing 6j-symbols from the diagram by encircling the vertices and 
crossings. This is useful in keeping track of the labels. Similarily one could also encircle the edges and faces 
to obtain the 3j and Ij-symbols. 

O'cro'' 

(r, a,b,c) = ^ JJ^^ 6(0, 1), r,)e(0, c, r,)e(li.c. rj)e(rj, b, r3)e(n , o, rj) 

Since the sum is over admissible colors only, we see that in this special case there is quite some simplifica- 
tion. The fact that the outside color is forces ri = b, r2 = c. Actually one can check that after substituting 
the formulas for the 6j, 3j and Ij'-symbols one gets that the answer is unless a = and 6 = c in which 
case we get a twisted unknot. 

2.4. Integrality. In [Cos09] it was shown that (r,7) € AfZ[A^'^] for some / depending on (r,7). In this 
subsection wc use the shadow formula to extend this result to prove that actually (r,7) e 
Here Qij) is a quadratic form depending on F in the labels 7, that takes values in Z/4Z. The proof below 
will indicate how to obtain an expression for Q(7) from an admissible coloring of the faces of a diagram of F. 
The form itself is independent of the choice of diagram, and coloring. By identifying Z/4Z with {0, 1, 2, 3} 
we have defined ((F,7)) = A^'^^'''' (F, 7) precisely. 

Proof, (of Proposition 1.1). Note that by Lemma 2.3 we can restrict ourselves to trivalent F. Pick any 
blackboard framed diagram D of (F, 7) and consider applying the shadow formula (15). Since by Costantino's 
result [Cos09] we already know the evaluation is a polynomial in A we can look at the degree of the terms 
in the sum. 

Now define (j) : Z[A=^^] Z/4Z to be a function such that P{a) £ ^'^(^)Z[A=^'']. As a first step we calculate 
(p for the building blocks of the shadow formula. By dividing the leading powers of the denominator in the 
formulas for the Gj-symbol we see that: 

(j){[k]) = 2k -2 
(j){e,a,b,c) = a^ + &^ + c^ 

i i<j 
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It suffices to sliow tliat every term in tlic shadow formula has the same value of (modulo 4). In order to 
checlc this we consider the effect on of a term when we increase one of the variables r by two. This is 
sufficient since a simple argument shows that any state can be reached from any other state by repeatedly 
increasing or decreasing one of the variables by two. 

So let 0(r) be the value of on a particular summand in which we ignore the terms not containing r. 
With respect to the region labeled r wc make a distinction between the edges and regions directly adjacent 
to it (notation: e\r or r^jr) and those edges and regions transverse to the region r (notation: e _L r). See 
also Figure 6. 



r2 



Figure 6. The region labeled r. The region ri the edge e are transverse to r. The region r2 and the edge 
e' are adjacent to r. 



Looking at the above formulae for </> and counting the crossings crossings next to r with sign a, we get 
the following expression for 4>{r). Here we abused the notation to make ri denote both the face and its value 
and (j(ri) is the sign of the corresponding crossing between r and r^. 

0(r) — r^(#adjacent edges + ^adjacent vertices + -^crossings + ^/^unsigncd crossings)+ 
r(2 + 2 ^ 7(e) + 2 + ^ 7(e) + J2 

Ti + cr(ri)) 

e|r 'fi\r e^^r ViA^r 

We will show that A = 0(r + 2) — 0(r) is divisible by 4. Since modulo 4 we have 
A = 2((r + l)#crossings + ^ 7(e) + ^ ''i + cf{ri)) mod 4 

we need to check that the term in brackets is even. To do so we use admissibility and write congruences 
modulo 2. By admissibility of the labels at the trivalent vertices and admissibility of the thetas coming from 
the edges in the shadow formula respectively we have: 

7(e) = 7(e') + 7(e") J'i = 7(^1) + 7(62) + r mod 2 

Here e, e', e" are the edges at a vertex and ei, 62, r, ri are the edges and two opposite regions at a crossing. 
Summing all these identities we see that every edge adjacent to r appears exactly twice, which shows that 

7(e) + ri = r^crossings mod 2 

Finally X^r _Lr '^(''i) = #crossings, therefore A is divisible by 4 and the proof is complete. □ 

3. Proof of Theorems 1.3, 1.4, 1.5 and 1.6 
3.1. Behavior of g-holonomic sequences under differentiation. In this section wc prove Theorem 1.4. 

Proof. Consider a q-holonomic sequence /„((?) G ^K?) that satisfies a linear recursion relation 

d 

^aj{q'\q)f^+j{q) =0 
j=0 
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where aj{u, v) £ Q[u, v] for j ~ 0, . . . , d. In the proof, the choice of coefficients Q[u, v] versus K[u, v] for a 
number field K does not play an important role, and therefore we assume K = Q. Differentiate with respect 
to q. We obtain that 

d d d 

(f6) 5^a,(q^g)/;+,.(g) +ng"-l^a,■„(g^g)/„+,(g) +5^a,,.(g",<7)/„+,(g) = 

i=o i=o j=o 

where aj^u = daj/du, aj^y — daj/dv and f'{q) ~ df{q)/dq. Recall now that the product and the Q[q]-linear 
combination of two g-holonomic sequences is q-holonomie; see [Zei90, WZ92]. Lemma 3.1 implies that for 
every J, the sequences {nq"~^aj,u{q" iq)fn+jiq)) and {aj,v(q"',q)fn+j{q)) are q-holonomic. It follows that 
the second and third sum in Equation (16) is (/-holonomic, i.e., it satisfies a linear recursion relation with 
coefficients in Q[q",q]. Substituting in this linear recursion the first sum of Equation (16), it follows that 
the sequence f^iq) is (/-holonomic. □ 

Lemma 3.1. For every fixed integer c, the sequence {nq'^") is g-holonomic. 

Proof. It is easy to verify that the sequence fn{q) — nq"^ satisfies linear recursion relation 

fn+2iq) - 2q'-'fn+iiq) + q^'fniq) = 

□ 

3.2. Evaluations of g-holonomic sequences at roots of unity are holonomic. Consider a q-holonomic 
sequence fn{q) G ^[9*^] that satisfies a linear recursion relation of the form 

(17) Crf(g", g)/„+d(g) + • • • + co(g", g)/„(g) = 

for all n where Cj{u, v) S Q[m, w] are polynomials in two variables for j = 0, . . . , d. The field Q of coefficients 
does not matter here and can be replaced by C without any change. Fix a complex root of unity The 
idea of the proof consists of a lucky ease, and a general case which reduces to a lucky ease after sufficient 
differentiation. The differentiation trick appears in an efficient algorithm to compute the Kashaev invariant 
of some knots in [GL05], and was also inspired by conversations of the first author with D. Zagier. Let us 
first discuss the proof of Theorem 1.5 when C = 1- In the lucky case, the set 

S^{j\0<j<d, c,(l,l)^0} 

is non-empty. Expand (17) as a power series in g — 1. The vanishing of the constant term implies that /«(!) 
satisfies a non-trivial linear recursion with constant coefficients 

^c,(l, !)/„+, (1) = 

Thus, the sequence (/n(l)) is holonomic. In the general case, there exists a unique natural number s > 
such that Cj{q"-, q) = ^j{q — lY + 0{{q — l)"'^^) and some 7j is nonzero. In other words, Cj{q",q) vanish to 
order s — 1 at g = 1 for all i, and some Cj{q^^, q) does not vanish to order s ai q ~ 1. In that case, observe 
that 

jes' 

where 

S' = {j \0<j< d, 7, ^ 0} ^ 

and c"f\i, f ) G Q[n\. This concludes the proof of Theorem 1.5 when C = 1- 

Suppose now that C is a complex root of unity of order TV. The problem is that cj^''(C":C) is no longer 
a polynomial of n even when s = 0. To overcome this obstacle, we consider n to be in a fixed arithmetic 
progression modulo N. In other words, fix i with < i < N ~ 1 and replace n by Nn -1- i in (17). Then, 
Cj{q^"~^^,q) can be expanded in powers of g — C with coefficients in Q(C)[7t.] for all i and j. In other 
words, for all i,j we have Cij{q) := Cj{q^"'^''' , q) G Q(C)[n][[(7 — C]]. Let us consider the column vector 
Xn = (/Arn(C), ■ • ■ jNn+N-iiO) G C^. If Cij(C) ^ 0, thcu fNn+i+j{C) is a Q(C) [n]-linear combination of 
fNn+i+j' iO for j' 7^ j. Of course, fNn+i+jiC) is the fc-th coordinate of Xn where k = i + j mod N. 
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Now we consider two cases. In the lucky Case 1, the foUowing set 

S = {{i,j)e[0,N-l] X [0,d]| c,,,(C)^0} 

is non-empty. S defines a graph G{S) defined as follows. It has vertex set {0, 1, . . . , iV — 1}, and an edge e 
between vertices k and k' iff there exist G S such that k = i + j mod N and k' = i + j' mod N. 

We consider two subcases. In the very lucky Case 1.1, G{S) is connected. The above discussion implies that 
the vector Xn and finitely many of its translates, put together in a column vector y„ satisfy a first order 
linear recursion of the form 

yn+i = A{n)yn 

for a square matrix A(n) with coefficients in the field Q(C)(?^). Lemma 3.2 below implies that every coordinate 
of y„ is holonomic. It follows that {fNn+i{C)) is holonomic (with respect to n) for every fixed i G [0,N —1]. 
Since (a„) is holonomic if and only if (aNn+i) is holonomic for all « G [0,7V — 1] (see [PWZ96]), it follows 
that fniC) is holonomic in the very lucky Case 1.1. 

In the not-so-lucky Case 1.2, the set S is non-empty but the graph G{S) is disconnected. In that case, 

differentiate (17) once and consider the pair of column vectors x'^ = (/^^(C), ■ • ■ >/wn+Ar-i(C)) S for 
r — 0,1. Consider also the set 

5'- = {(*,j)e[0,iV-i]x[o,d]|cg(C)^0} 

for r = 0, 1. S'^US^ gives rise to a graph G'(5"u5^) with vertices pairs (r, k) for r = 0, 1 and fc = 0, . . . , A^— 1. 
There is an edge between (r, fc) and (r', fc') if there exist £ , G such that k = i + j mod 

and k' = i + j' mod N. Note that by assumption, S*" = 5 is non-empty, and there is a natural inclusion 
of G{S) into G(5° U S^). In the not-so-lucky Case 1.2, the graph G{S^ U S^) is connected. The argument 
of Case 1.1 shows that fn{C) is holonomic. If G(5'° U S^) is disconnected, differentiate (17) once more. It 
is easy to see that after a finite number t of differentiations, the graph G{S'^ U . . . 5*) will be connected 
and consequently fn{C) is holonomic. Differentiating (17) t times and using induction implies that /n(C) is 
holonomic. 

This concludes the proof of Theorem 1.5 when S is non-empty. In the unlucky Case 2 where all functions 
Cij{q) vanish to first order ai q ~ (, there is a natural number s such that cfj\C) = for s' < s and all 
i,j and in addition c\'^j{C) 7^ for some In that case, replace S by the set 

S' ^{ii,j)\0<j<d, 0<i<N- 1, cg(C) ^ 0} 

and repeat the above proof. 

Lemma 3.2. Consider a sequence of vectors x„ E C' that satisfy a first order linear recursion of the form 

Xn+i = A{n)x„ 

where A{n) is an ^ x ^ matrix with coefficients in C[n\. Then every coordinate of Xn is a holonomic sequence. 

Proof. Consider the operators u and N which act on a sequence (a„) by 

(jya)„ = na„, {Na)n = a„+i. 

These operators satisfy the commutation relation Niy = lyN + N and generate an Ore ring R = <C[i']{N). 
The matrix A{n) defines an ii- module M as in [vdPS03]. The lemma follows from the cyclic vector theorem 
applied to the module M; see [vdPSOS, Prop. 2. 9] and also [Kat87]. The proof is constructive and can be 
implemented for example in [AB06]. □ 

This finishes the proof of Theorem 1.5. □ 
Note that Theorem 1.5 fails when q is not a root of unity. For example, consider the g-holonomic sequence 

o-nil) = and fix g = cj, a complex number which is not a root of unity. Suppose that the sequence 6„ = 

is holonomic. So for all natural numbers n we have 

d 

^Cfc(n)6„+fe = 

k=0 
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where Cj{n) G Q[n\. Dividing by 6„ we get 

k=0 

Collecting the coefficients of a fixed power of n, we find that there are Ci . . . Cd G C such that 

D 

J2 Ck^'^'' = 

fe=0 

In other words, for all n, is a root of the polynomial X^feLo ^kX^ ■ It follows that = 1 for some 



m. 



3.3. Proof of Theorem 1.6. In this section we prove Theorem 1.6. Recall from Remark 1.7 that this 
theorem is special to the quantum spin network evaluations and not valid for general (/-holonomic sequences. 
The main ideas that go into the proof of Theorem 1.6 are the following: 

(a) its validity for the local building block of quantum spin networks (namely the quantum 6j-symbol), 

(b) The shadow formula (15) for the evaluation of an arbitrary quantum spin network where the sum- 
mand is roughly a product of quantum 6j-symbols divided by quantum 3j-symbols, 

(c) a lemma for the growth rate of the coefficients of inverse quantum factorials, in the spirit of [GL05, 
Thm.15]. 

Let us begin with the first ingredient. 

Lemma 3.3. Theorem 1.6 holds for the quantum 6j-symbols. 

Proof. Equation (12) shows that /«((?) = {/%.,n'y) is a Laurent polynomial in q^/^ which is written as a 
1-dimensional sum. Let ||/i(q)||i denote the sum of the absolute values of a Laurent polynomial in g^/^. It 
is well known that the quantum binomial coefficient (10) is a Laurent polynomial in q^f^ with nonnegative 
integer coefficients. Moreover, its evaluation at g = 1 is the usual multinomial coefficient. The sum of the 
latter is at most r°. It follows that 

a! 



fli! . . . a^! 



_ai, 02, . . . , ar_ 

Equation (12) and the above inequality conclude the proof of the lemma. □ 

Now consider an arbitrary quantum spin network (r,7) and its evaluation /Nil) — N^){q). Viewing 
Ij-symbols as special cases of 6j-symbols and crossings as a power of q times a 6j-symbol, the shadow 
state-sum formula (15) expresses fNil) by 

M{k) 



(18) fNiq)= J2 



Dik) 

where the summation is over the set of lattice points of NP (for some rational convex polytope P) and M{k) 
is a product a power of q and of quantum 6j-symbols at linear forms of k, and D{k) is a product of quantum 
factorials at linear forms of k. The fact that D is a product of quantum factorials only follows from the 
formula for the 3j-symbol, see Equation (11). We know that /7v(<z) & Z[q^^^^] by Proposition 1.1. Theorem 
1.5 implies that the exponents of q in (F, Nj){q) are bounded above and below by quadratic functions of N. 
Let {q)n ~ nj=i(l ^ 9"') for n e NU {oo}. It is well-known that 

^ oo 
'''^>°^ k=0 

where € N is the number of partitions of k; see [And98]. Moreover, for every k and n, the coefficient of 
g*^ in l/{q)n is a natural number which is at most pk- In addition, we have 
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where C = 7r-y/2/3; see [And98]. It follows that the coefficient of in 1/D{k) is exponentially bounded by 
n^, and since n = 0{N'^), the coefficient is exponentially bounded by N. Since the number of fc-summation 
terms in the state-sum (18) is bounded by a polynomial function of TV, this completes the proof of Theorem 
1.6. □ 

Remark 3.4. The above proof is similar in spirit with the proof of [GL05, Thm.l5]. 

3.4. Proof of Theorem 1.3. Fix a quantum spin network (r,7). It follows from the shadow state sum in 
Section 2.3 and Lemma 2.3 that the standard evaluation (F, n'y){q) is a sum of a balanced g-hypergeometric 
term. Hence, by the fundamental theorem of WZ theory (see [WZ92]), it follows that {r,n'^){q) is q- 
holonomic. Since {{T,nj)){q) = q'^^'"'^^^^{T,n'y){q) and Q{n'y) is periodic, ([T^wy^lq) is g-holonomic as well. 

Fix a complex root of unity C,, and a natural number r > 0, and let a„ = ((F, n7))(''^(C). We will apply 
Theorem 1.5 to check conditions (a)-(d) of Definition 1.3. 

(a) a„ € Z[^] as follows from Proposition 1.1. 

(b) (b) Theorem 1.6 shows that the coefficients ((F, ?^7))^f )(C) of the generating function are exponentially 
bounded. The same goes for their Galois conjugates since the action of the Galois group can only 
replace the root of unity ^ by another root of unity, so Theorem 1.6 still applies. 

(c) Since an € are algebraic integers, their common denominator is 1. 

(d) By theorems 1.4 and 1.5 the coefficients of the power series are holonomic, hence the series itself is 
holonomic as well. 

This concludes the proof of Theorem 1.3. □ 



4. Examples 



In this section we consider various simple examples to illustrate the theorems. For the sake of simplicity 
we surpress the function A~^'^'^^ so we choose to work with (F,7), rather than the more correct ((F,7)). 

4.1. Evaluations of the q- multinomial coefficients at a root of unity. 
Lemma 4.1. Let q = and set ai = AiN + ai. We have 



(19) 



ai + . 

01,02, 



+ an 

■ • 1 



«! + ...+«„ 
"1, "2, ■■-,«„ 



id) 



This was proven for Gaussian (asymmetric) binomial coefficients in [Des82]. The extension to the multi- 
nomial case is straightforward and the power of q appears when one converts to symmetric multinomial 
coefficients. Note that the power of q is actually a sign since aja^ — ajUk ~ —N{AjAkN + ajA^ + Ajak). 

Now consider the growth rates of the theta evaluations. Fix and admissible triple a,b,c E N, and set 
^ ^ -a+b+c ^ Q ^ CL^^ C = ii+l^ then we have 



(9, an, bn, cn) = ( ( i, n{A + B + C)) 

At g = 1 we find using Stirling's formula (sec [01v97]) 

,A+B+C (A + B + C) 



n{A + B + C) 
nA, nB, nC 



Ai 



A+B+C 



i-iy 



At the A^-th root of unity we find 

(0, an, bn, cn) {e~ ) ~ Pin) 



AAbBcC 



I I I I I I I "C I 

L AT J L AT J L AT J 

I ""^ I I I I nC I 

L w J ' L AT J 7 L AT J 



where P{n) is a periodic function in n of period 2A^. Another application of Stirling shows that Ajv C 
f^2Ar|Ai| « , where fiw denotes the set of all A-th roots of unity. 
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4.2. The g-diflference equation of the regular 6j-symbol. Wc compute a second order recursion relation 
for the regular quantum 6j-symbol 

an{q) = {/:, ,ni){q), 7 - (2, 2, 2, 2, 2, 2). 

using the g-WZ method, implemented in Mathematica by [PR] and used as in [GvdV09]. The recursion has 
the form 

2 

(20) ^Cfc(q,g"K_fc(<z)-0 

k=0 

where 

Co(g, q") = {q - (-1 + g")'' {q + q") {4q^° + Aq^"" + 6g^+" + 4g«+'" - g«+'" + 5q' 

_6g5+4" - 2g'^+'*" - 65^+"*" - 4g''+^" + 2g^+'" - 4g'^+^" - 65^+^^" - 2q*+'''" - 6q'' 
+5g^+"" - q^+^" + 4q^+'" + 69^+^") 



/ n\ —1 + 671 / n\3 / 2ri\ / a 14 ^ 15 , 18n , ^ 19n , t r\ 20n , in 21n , r, 22n 

ci(q,g ) = -q {q-q) [q~q ) (-4g - 4g + g + 7g + lOg + lOg + 7g 

+g2^" - 10g"+" - 8g"+" - 10g^''+" - 10g^'+'" + 6g"+'" + 6g"+'" - 10g^''+'" + q^°+^" 
_9^ii+3n _^ ^i2+3n ^ 28g"+'" + g"+'" - 9q'''+''" + gi''+'" + 7g''+''" + 4q'°+^" + 14g"+''" 

+36g"+^" + 36g'^+'^" + 22g"+^" + 2gi^+^" + 10g'®+^" + 10g''+''" - 3g*'+'^" - 16g^+''" + 4gi°+''" 

-I n 11 + Gn 1 n 12+6n . a 13+6ri 1 14+6n o 15 + 6n . t ri 16 + 6n . ^ 6 + 7n o 7+7n . 8+7n 

— 17g — 17g + 4g — Ibg — 3g + lOg + 7g — 3g +g 

, a 9+7n o 10+771 . ll+7n o 12+7n , ^ 13+7n , 14+7n r, 15+7n , ™ 16+7n , 5+8n 

+Dg — 3g +g — 3g +Dg +g — 3g +7g +g 

_-j^gg6+8n _ 3]^g7+8n _ .^^8+8n _ 4^^9+8" _ 37^^"+*" „ 37g"'""'"^*" — 47g"'"^^®" — r,^13+8n _ 3^^^14+871 
_^ggl5+8. ^ gl6+8. _ g^5+9n _ ^g^6+9n ^ ^^7+9n _ ^g^8+9n _ ^2g9+9n _ g^lO+g. _ ^^q^Wn 
_18ql2+9n ^ 4gl3+9n _ I8gl4+9n _ g^l5 + 9n _ g^4+10n _ g^5 + 10„ ^ g^6+10n _ ^g7+10n ^ 13g8+10n 
+ 13g®"*'^°" + 13g^'''^^°" + 13g"^^"''^''" - y^l2 + 10n _|_ g^l3 + 10n _ g^l4+10n _ g^l5 + 107i _ 4^3 + 1171 

-2g''+"" - g^+"" - 13g'+"" - 24g'+"" + 2g«+"" - 43g^+"" + 2g^''+"" - 24g"+"" - 13g''+"" 
_qi3+iin _ 2^14+1177 _ 4^15+1177 ^ iig3+i2n ^ 31g^+^^" + 21g^+^^" + 449=^+^^" + 53g^+^^" 

+56g«+^^" + 56g^+^^" + 535^"+^^" + 44g"+^^" + 21g^^+^^" + 31g"+^^" + llg"+^^" - 4g^+"" 
+ 18g^+'^" - 12g^+"" - 14g''+"" - 26g''+"" - 40g^+"" - 36g**+"" - 40g^+"" - 26g'°+'='" 
_14q"+i3" _ I2gi2+i377 ^ 18^13+1377 _ 4^14+1377 _^ 1^^2+1477 _^ 31^3+1477 _^ 21g4+i^" + 44g^+"" 

+53g«+"" + 56g^+"" + 56g«+"" + 53g'-'+"" + 44g^°+^'*" + 21g"+"" + 31g^^+"" + llg^^+"" 
_4gi+i5" _ 2g2+i5" - g^+15" - 13g^+'''" - 24g"+''" + 2g''+''" - 43g"+''^" + 2g*+''" - 24g'+''" 

— 13g"'"'''''"'"^" — _ 2q^^'^^^" — 4^"'"''+^^" — gg"""^^®" — Sg'^^"'*'" + Qg^^"*""" — 7^*+-'-®" 13g^^"'"^" 
+13g^+^«" + 13g"+^''" + 13g**+^«" ~ 7g^+^«" + 9g^°+i«" - 8g"+^'^" - 9gi^+^«" - 8g^+"" - 18g^+i"" 
+4g3+i7n _ 18^4+1777 „ i2g^+^^" - 8g''+^"" - 12g^+^^" - 18g«+^"" + 4q''+"" - 18g"+"" - 8g"+^"" 
-18g^"''^*" — 31g^"*'^*" — 7^^+18" _ 47g4+i8" _ S7q^^^^" - 37g'''^"^*" - 47g^''"^**" - 7^8+18" _ 3-1^^9+1871 

10 10+18n , ll+18n o l+19n , 2+1977 . e 3+1977 o 4+19n , 5+19n q 6+1977 . ^ 7+1977 . 8+1977 

-18g +g — 3g +g +6g^ -3g +g -3g^ +6g^ +g 

_3g9+19n + 7gl0+1977 _ 3^1+2077 _ ig^2+2077 4^3+2077 _ ^^^4+2077 _ ^^5+2077 _|_ 4^6+2077 _ ^Q^7+20n 

_3^8+20n + 10g«+2"" + 2g^+2^" + 22g2+2^" + 3Gq'+^^" + 36g'+'i" + 22g''+2^" + 2g'+'^" + 10g'+'^" 

, . l+22n , 1 . 2+2277 , , o 3+2277 , , o 4+2277 , , . 5+2277 , . 6+2277 . - 7+2277 ri 1+2377 . 2+23n 

+4g + 14g + 13g + 13g + 14g +4g +7g — 9g +g 

_^28q3+23n ^ g4+23n _ gg5+2377 ^ ^6+2377 _ + Qg^+^-i" + Qq'+^^^ - 10g*+"''" - 10g'+'"" 

_8g2+25n _ 10g3+25- _ 4^1+2677 _ 4g2+2677N 
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/ n\ 1 /t I 7i\ / 2 3n\'i / 4 3n\*^ / 2 . 2n , l + / . 377, „ 4n . 577, „ 5n 777, 

c2(g,Q ) = -^{i + q )[q -q ) [q -q ) [q +q +q ) (-4g + g +6q +4q +6g +q 

-6g'+" - 4q'+2" - 5g'+^" + + 2g'+^" + 6^''+''" - 2g'+^" + 4q^+'" + 2q'+''" + 6g'+^" 

-5g^+^" + - 4q'+'" ~ 6q'+^" - 4g^+^"") 

The above recursion can be used to deduce a recursion relation for the sequence a„(C) at every fixed root 
of unity following and illustrating the proof of Theorem 1.5. For example, to obtain the recursion of the 
sequence a„(l), use: 

co(g) = -4((-l + n)3n5(106 - 230n+115n2))(g-l)i" + O((g-l)") 
ci(g) = 4(-l + nfi-l + 2n)(-864 + 4470n + 211471^ - 5100371^ + 12008971"^ - 11350571^ 
+37835776)(g-l)i° + 0((g-l)") 

C2{q) = -3 24((8 - 1877 + 9772)4(-9+115n2))(g- 1)10 + 0((g- 1)11) 

An alternative method to obtain the recursion relation of the sequence a„(C), using the WZ method, is 
described in the next section. 



4.3. The regular 6j-symbol at roots of unity. The regular 6j-symbol corresponds to the labeling 7 
(2,2,2,2,2,2) and set 

a„ = (^4,77 7)(1), 6„ = ( ,;K,n7)(-l). 

Equation (12) gives 



k—3n 



(fc- 377)14(477 -/c)! 



Using the WZ method (implemented in Mathematica by [PR] and using it as in [GvdV09]) it follows that 
(a„) satisfies the recursion relation 

-81 (2 + 3n)' (4 + 3n)'' (451 + 460n+115n^) a[n] + 

(1+n)^ {3 + 2n) (319212 + 1427658n + 2578232 n^ + 2423109 n^ + 1255 139 n" + 340515 + 37835 n') a[l+n] - 
(1 + n)^ (2 + n) = (l06 + 230n+115n2)a[2+n] ==0 

This linear recursion has two complex conjugate formal power series solutions a±^n where 

1 / -304- 3n\/2 25879 + 18352^^/2 71176912 + 632307HV2 
- + 57677 + 33177677^ + 573308928773 

5(-2742864803 - 10265264480^^^) -82823449457840 - 12750219708659i\^ 
660451885056774 38042028579225677^ 
, -61273664686901989 + 213495822835779152iV2 1 



657366253849018368776 



77' 



When q = —1. Equation (12) and Lemma 4.1 imply that bn — for odd 77, and 

4n 

hn = E (-1) 



fc! 

(fc- 377)14(477 -fc)!3- 



It follows that the sequence (6271) satisfies the following recursion relation 

81 (l + 3n)' (2 + 3n)' ( 3822 + 10 7 7 9 n + 11 352 + 52 90 + 920 n' ) a[n] - 

1 (1 + n)^ (41796 + 522312n + 27 92 052 n2 + 8292 059 n^ + 15002 105 n' + 171587 65 + 12460766 n' + 5571784 n' +13998 95 n' + 151340 n' 
a[l+n] + (1+n)^ (2+n)^ (25 + 265n + 1002 n^ +1610 + 920 n') a[2 + n] == 
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This linear recursion has two complex conjugate formal power series solutions b± 2n where 



1 



■(329 + 460i\/2)' 



1 



-472- 19i\/2 105199 + 20632^72 22386472- 9304565^72 



^5/2 • — , ^ , g^g^ , 33i776n2 ' 573308928n3 

-30711007135 - 48640734448iv^ -80744339543960 + 2822061829369iv^ 



660451885056n4 
17678244315725891 + 219394408835134568iv^ 



380420285792256^5 



657366253849018368n 



6 



o 



1 



5. Open problems 



Fix a cubic ribbon graph F with edge set i?(F), a complex root of unity C, and consider the generating 
sevies 

^r,c=E«r,7))(C) n &m)[[^e.eeE{T)]] 

7 ee_E(r) 

In [GvdV09, Thm.5], we proved that when C = 1, St,i is a rational function, i.e., belongs to the field 
Q(ze,e € EiV)). It was explained in [GvdV09, Thm.6] that the generating series -Fr,7,c (4) is a diagonal of 
St,q, and consequently the rationality of St.q implies that is a G-function. The rationality of St,q 

follows from the so-called chromatic evaluation of classical spin networks. For details, see [GvdV09] and also 
[Wes98] . The chromatic evaluation seems to break down in the case of complex roots of unity. 

Question 5.1. Is it true that for all cubic ribbon graphs F and all complex roots of unity, is a rational 
function, i.e., belongs to the field Q{C){ze,e G ^^(F))? 
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